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Solving the one-dimensional time-dependent Schrodinger 
equation for simple model potentials, we investigate 
resonance-enhanced high-order harmonic generation, with 
emphasis on the physical mechanism of the enhancement. By 
truncating a long-range potential, we investigate the signif- 
icance of the long-range tail, the Rydberg series, and the 
existence of highly excited states for the enhancements in 
question. We conclude that the channel closings typical of 
a short-range or zero-range potential are capable of generat- 
ing essentially the same effects. 
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I. INTRODUCTION 

A very intriguing feature of above-threshold ionization 
or high-order harmonic generation concerns the depen- 
dence of these phenomena on the intensity of the driv- 
ing field. The photo-electron or high-order harmonic 
peaks, as functions of the driving-field intensity, present 
resonance-like enhancements, such that a variation of a 
few percent in the external-field strength may drive up 
the spectral intensity by an order of magnitude. These 
enhancements have been observed experimentally by sev- 
eral groups for above-threshold ionization (ATI) 
and, recently, for high-order harmonic generation (HHG) 

. A concomitant effect in ATI is a variation of the con- 
trast of the spectrum [Q. 

Early numerical observations that enhancements of 
ATI go hand in hand with enhancements of HHG were re- 
ported in Ref. |^ . The existence of enhancements for both 
phenomena is not surprising, and is related to their com- 
mon physical origin. Indeed, HHG and ATI present very 
similar spectral features, which are explained by similar 
physical pictures. These features are a wide energy range 
with approximately equally strong harmonics or photo- 
electron peaks, known as "the plateau", followed by a 
sharp decrease in the harmonic or photo-electron signal, 
known as "the cutoff". HHG is described by the so- 
called "three-step model" , in which an electron is ionized 
through tunneling or multiphoton ionization, is acceler- 
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ated by the field and driven back to its parent ion, where 
it recombines to the ground state, emitting its energy as 
one harmonic photon [^j^. A similar process is respon- 
sible for the plateau in above-threshold ionization, with 
the main difference that, instead of recombining with the 
parent ion, the electron is elastically rescattered off it [|| . 
In either case, the precise shape of the atomic potential 
is not very important, and it can be approximated, for 
instance, by a zero-range potential. This approximation 
describes very well the spectral features near the high- 
energy end of the plateau. 

The resonance-like enhancements, however, primarily 
occur in the first half of the plateau, for which both the 
external driving field and the atomic potential are ex- 
pected to influence the harmonic or photo-electron emis- 
sion In fact, very different arguments have been put 
forward to explain these enhancements. Several studies 
attribute these features to Rydberg states that, for an 
appropriate ponderomotive upshift, become multiphoton 
resonant with the ground state. Namely, a free electron 
in a laser field acquires a field-dependent energy shift 
by the ponderomotive potential Up — (>l^(t))j /2m, 
where A{t) is the vector potential of the laser field and 
denotes the average over a field cycle. Highly ex- 
cited Rydberg states tend to undergo about the same 
shift as free electrons [0. The result of such a multi- 
photon resonance is either an increase in ionization, or 
the electronic wave packet is trapped near its parent ion 
for relatively long times, originating resonance- like struc- 
tures in the spectra [ pH . The mechanism is quite similar 
to the Freeman resonances , which dominate the low- 
energy ATI spectrum where rescattering plays no role. 
According to this physical picture, the presence of high- 
lying Rydberg states is essential for the existence of the 
enhancements. 

An - at least at first glance - completely different view 
relates these enhancements to channel closings that, by 
the same ponderomotive-upshift mechanism, may move 
into multiphoton resonance with the ground state. If, 
due to this shift, N photons are no longer sufficient for 
the electron to reach the continuum, one refers to the 
A''-photon channel as having closed. At an intensity cor- 
responding to a channel closing, the electron is released 
in the continuum with a vanishing drift momentum. In 
consequence, in the course of its oscillatory motion in 
the laser field, it will return many times to its parent ion 
and upon each revisit have the opportunity to rescat- 
ter. Quantum mechanically, the corresponding probabil- 
ity amplitudes interfere, and a constructive interference 
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manifests itself as an enhancement 0. Such an effect 
does not require the existence of excited states or Ryd- 
berg states so that the atom can be modeled by a binding 
potential of zero range, which neither supports excited 
bound states nor resonances in the continuum P, p^ , p^ . 
The zero-range potential affords the concept of "quan- 
tum orbits" which allows for an almost analytical ap- 
proach to intense-laser-atom phenomena [p^ . However, 
one has to keep in mind that real atoms do have long- 
range-potential tails, so that channel closings are diffuse 
owing to the presence of the Rydberg series. In addition, 
they support various bound states, whose influence on 
the enhancements is not entirely clear. 

In this paper, we perform a systematic study of the 
influence of both the laser field and the atomic potential 
on these enhancements, for high-order harmonic genera- 
tion, by means of simplified atomic models. We address 
the question of which one of the existing physical inter- 
pretations is ultimately correct or whether both pictures 
are complementary aspects of a more complete descrip- 
tion. In particular, we investigate the importance of the 
highly excited states in the process, and whether or not 
they are crucial for the feature in question. 



II. TRUNCATED SOFT-CORE POTENTIALS 

We compute the harmonic spectra using the numerical 
solution of the time-dependent Schrodinger equation 



^i\m) 



+ V{x) -p-A{t) 



for a one-dimensional model atom |lq| initially in the 
ground state of a binding potential V{x) and subject to 
a laser pulse with the field E(t) — —dA{t)/dt. Atomic 
units are used throughout. We consider a monochromatic 
laser field 



E{t) = Eq sinwi. 



(2) 



and the harmonic spectra are calculated from the dipole 
acceleration i = {i}j\t)\-dV{x)/dx + E{t) \i){t)) 0- We 
take the smoothly truncated soft-core potential 



V{x) = 



with 



fix) 



1 (|x| < ao), 

(ao < \x\ < L), 



(3) 



2(L-ao 

(|x| > L), 



(4) 



so that V{x) = for |a;| > L. We choose oq of the order 
of a few atomic units, and L of the order of the electron 
excursion amplitude, so that L = ra = tEq juP' with the 
parameter r of order unity. Setting f(x) = I in (^) gives 



the untruncated soft-core potential. By an appropriate 
choice of the parameters L and ao, it is possible to alter 
the highly excited bound states leaving the ground state 
and the low excited states practically unaffected. 

Let us assume that an atom, initially in the ground 
state with energy Eq, is ionized by iV photons of frequency 
cj, such that the electron reaches the continuum with the 
lowest energy possible, that is, with a drift momentum 
(outside the range of the binding potential) of zero. The 
energy of the N photons must account for the binding 
energy and the kinetic energy of the oscillatory motion 
(the ponderomotive energy [/p), so that 



(5) 



For intensities slightly larger than specified by the condi- 
tion (H), at least TV -I- 1 photons will be necessary for ion- 
ization, such that Eq. (^) defines the A^-photon channel- 
closing intensity. The intensities that solve the channel- 
closing condition (||) form a comb whose teeth as a func- 
tion of 7^ = Up/uj are separated by unity. If there is an 
excited bound state with the (field-free) energy e„ and 
if this state undergoes the same ponderomotive upshift 
as the continuum [ pO[ , then multiphoton resonance with 
the ground state occurs for intensities such that Eq. (|^) 
is satisfied with replaced by Eq — Sn, 



Up = Nuj. 



(6) 



For a long-range potential [such as, in one dimension, 
our untruncated potential (^)] the true continuum is pre- 
ceded by the Rydberg series so that one may question 
the significance of the channel-closing condition (^) for 
any physical phenomenon. For a finite-range potential 
[such as the truncated potential (^)] the Rydberg series 
is replaced by a finite sequence of bound states whose 
number decreases with decreasing L. 

Below, we will seek to answer the following questions: 
Is the very existence of an enhancement contingent on 
the shape of the binding potential? If an enhancement 
exists, does it occur at a channel-closing intensity (|^) for 
some N , or is related to a multiphoton resonance (§) with 
a certain excited state or is it unrelated to either? 
For a truncated potential, do enhancements occur both 
at the channel-closing intensities and at intensities where 
an excited state becomes multiphoton resonant? To what 
extent does the harmonic spectrum depend on the shape 
of the potential? 

For all examples presented in this paper, we choose 
/? = 2.1 and a — 0.2 in the soft-core potential (|^). In 
this case, the field-free energies of the ground state and 
the first four excited states are listed in Table || for the 
untruncated potential as well as for various truncations. 
The Table shows that the ground-state energy is virtu- 
ally unaltered by the truncations we consider while the 
excited states are more and more affected. For the un- 
truncated potential, the first excited states are followed 
by the Rydberg series. For all truncations considered, 
the states listed are the only ones that survived. By and 



2 



large, for the parameter range used, the truncation elimi- 
nates the Rydberg series, changes the energy of the third 
excited state, and leaves the more deeply bound states 
unchanged. The shape of the effective potential barrier 
Vctf — V{x)—xE{t) also remains very similar for all cases. 



III. RESULTS 

As a first step, we investigate a harmonic spectrum 
as a function of the laser intensity with regard to the 
existence of enhancements for particular intensities. The 
frequency of the laser field is taken as w = 0.07600 a.u., 
such that Eq. (||) predicts channel-closing intensities for 
near-integer values of 77 = Up/ui. The other intensities 
that correspond to resonances with excited states are also 
listed in Table |. For the frequency and the intensity 
range considered in this paper, the Keldysh parameter 
7 = ■y/|eo|/2t/p lies within the interval 1.66 < 7 < 0.86, 
which is mostly within the multiphoton regime. 
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FIG. 1. Harmonic spectra for ui = 0.076 a.u., Eq = 0.7566 
a.u., and several values of ?; = Up/ui, for the untruncated (a) 
and truncated (b) soft-core potential, with L ~ 31.78 a.u. 
and oo = 3 a.u. The parameter L was chosen as twice as the 
electron excursion amplitude a for rj — 4.8. The harmonic 
intensities are connected by lines to guide the eye. The thick 
lines correspond to the field intensity for which the enhance- 
ment is maximal. 

Figure |^ displays two harmonic spectra, computed for 
the untruncated (a) and a truncated (b) soft-core poten- 
tial. Both spectra are very similar, with pronounced en- 
hancements near the 13th harmonic. For the untruncated 
and truncated cases, we observe maximal enhancements 
at 7] = 4.3 and 77 = 4.1, respectively. These values do 
not coincide with the channel-closing intensity 77 — 4.04 
(modulo any integer) predicted by Eq. (|5|). This shows 
a clear influence of the binding potential on the abso- 
lute intensity for which these features occur. However, 
comparing parts (a) and (b) of Fig. |l| we notice a very 
remarkable fact: the spectra of the untruncated (a) and 
the truncated (b) potential are almost identical provided 
we compare a spectrum at the intensity r/ for the former 
with a spectrum at the intensity ry — 0.2 for the latter. 
This holds for all intensities considered. Table || shows 
that the truncated potential no longer supports the Ry- 
dberg series and its excited state |3) is very close to the 
continuum. We conclude that the Rydberg series has 
no visible impact on the shape of the harmonic spec- 
trum. Checking again Table | we further conclude that 
the intensities where the enhancements are maximal are 
compatible with Eq. with n = 3. 

Further support for these conclusions comes from Fig. 
H, where the yield of the 13th harmonic is plotted as a 
function of the scaled intensity 77 = Up/to. In this fig- 
ure, we also investigate the effect of the truncation in 
more detail, for a wider range of intensities and trun- 
cating parameters. The most striking feature is that the 
main effect of the various truncations is a rigid horizontal 
shift of the yield- versus-intensity curve. Remarkably, this 
statement includes the pronounced dips which are due to 
quantum interference between different quantum orbits 
0,[5|. While this may be plausible for the large values of 
the truncation parameter (L = 2a) where the truncation 
mostly affects the potential outside of the classical elec- 
tronic excursion, it is quite surprising for the small value 
L = 0.3a. A similar shift pattern is observed for all low 
plateau harmonics (not shown). 
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FIG. 2. Intensity of the 13th harmonic as function of 
rj — Up/iu. Part (a): Comparison between the untruncated 
and various truncated soft-core potentials; the two values 
of the parameter L, L = 31.78 a.u. and L = 4.77 a.u., 
correspond to 2a{ri = 4.8) and 0.3q(j7 — 4.8), respectively. 
The binding energies for the various truncated potentials are 
given in table ^ Part (b): Enlargement of part (a), for 
3.4 < rj < 5.4. The values of this parameter, for which the 
enhancements are maximal, are marked by arrows. 

Next, we investigate the quantitative amount of the 
afore- mentioned shift. Comparison with Table || shows 
that all of the enhancement peaks present in Fig. are 
compatible with the intensities predicted by Eq. (Q) for 
n = 3. As long as the state |3) exists as a bound state, 
we conclude that the enhancement is due to a multipho- 
ton resonance with this state, upshifted by the whole 
amount of the ponderomotive energy. However, the en- 
hancement still exists (near integer rj), though quanti- 
tatively smaller, in the case where the truncation (for 
L = 0.3a) has eliminated the third excited state. Here 
we have encountered a case of a pure channel-closing en- 
hancement which is due to a multiphoton resonance with 
the ponderomotively upshifted continuum threshold. 

Why does a multiphoton resonance with the third ex- 



cited state (as long as it exists) lead to a pronounced res- 
onance in the spectrum while a multiphoton resonance 
with other excited states does not? We can answer this 
question by inspecting the wave functions {x\n) of the ex- 
cited states. Please, notice that our calculations are still 
within the multiphoton regime of ionization. In this case, 
one expects that multiphoton resonance with an excited 
bound state is particularly relevant if the wave function 
of the latter is concentrated near the turning points of 
the wiggling motion of a classical electron with a drift 
momentum of zero, i.e., when l(a;|n)p has its maxima 
near « a = Eq/uP'. Figure || shows that indeed this 
is well satisfied by the state |3). 




X (a.u.) 

FIG. 3. The probability distributions |(a;|n)|^ for the first 
four excited states of the untruncated potential The 
range of excursion amplitudes a — Eo/lo'^ — 2-^77/0; that 
is covered in Fig. ^ (for 2.3 < rj < 6.3) is shaded. Truncation 
has a minor effect for n = 1 and 2. For n = 3 this effect is 
becoming noticeable, and for n = 4 it is strong. All higher 
excited states are eliminated by the truncations considered. 

It is not surprising that ionization is affected by the 
resonances in a similar fashion as HHG. In Fig. ^ we plot 
the normalization of the time-dependent wave function 
at the end of the pulse, as a function of -q for the trun- 
cation parameters used in Fig. |^. This normalization is 
smaller than unity, since part of the wave function is ab- 
sorbed by a mask function, used in our computations to 
eliminate spurious reflection effects. It is a good mea- 
sure of irreversible ionization, since the mask function is 
located at about ten times the electron excursion ampli- 
tude. In the figure, there are clear dips slightly preceding 
the intensities, for which the channel closings occur. The 
particular case, for which |3) is absent, agrees with the 
results in for a zero-range potential. 
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FIG. 4. Normalization of the time-dependent wave func- 
tion at the end of the pulse as a function of = Up/u), for the 
truncation parameters of Table |. The pulse length is eight 
cycles plus a linear turn-on of two cycles. 

There are still features in the enhancements, such as 
their strength and harmonic range, which depend on the 
atomic binding potential and are not yet completely un- 
derstood. In fact, studies performed by us for short- 
and long-range potentials of various shapes indicate that, 
even though the enhancements are a generic feature be- 
ing present for all cases, their strength, regularity and 
harmonic range depend on the potential in question. An 
interesting feature is the fact that the enhancements are 
weaker for short-range potentials. Since the tail of the 
potential plays no role in this respect, the effective po- 
tential barrier may influence these features. 

IV. DISCUSSION AND CONCLUSIONS 

The picture that emerges from these calculations is 
this: high-order harmonic generation can be strongly en- 
hanced by a multiphoton resonance with a ponderomo- 
tively upshifted excited bound state close to the contin- 
uum threshold. In the cases that we investigated the 
relevant (field- free) bound state \n) was the state where 
|2;n,max| ~ Oi = Eq / uj^ with Xn,m.ax the valuc of X that 
renders |(a;|n)p maximal. Bound states higher than this 
do not lead to noticeable enhancements, nor do the exis- 
tence or nonexistence of a Rydberg series have an effect, 
neither on the enhancements nor on the entire harmonic 
spectrum. In all cases considered, the relevant bound 
state was \n = 3). When this state was eliminated by 
sufficient truncation of the potential, the role of the reso- 
nant bound state was taken over by the ponderomotively 
upshifted continuum threshold. Hence, in all cases, the 
resonant intensities are well described by Eq. (|^), pro- 
vided we employ an effective continuum threshold that 



is given hy e = Eq — Sn where £„ is the energy of the cru- 
cial bound state if such a state exists and zero otherwise. 
Remarkably, if the value of ea changes owing to a change 
of the potential, the effect on the harmonic spectrum 
is largely a horizontal shift of the yield-versus-intensity 
curve of the various harmonics. 

This picture allows one, for the purposes of HHG, to 
model the atom by a short-range or zero-range potential 
using a binding energy that is adapted to the energy dif- 
ference between the ground state and the relevant excited 
state of the real atom that is supposed to be modeled. 
In this case, excited states need not be considered and 
the enhancement can be attributed to a multiphoton res- 
onance with the continuum threshold and the generation 
- in the three-step model - of an electron with a vanishing 
drift momentum. The picture of an effective continuum 
threshold is also intuitively appealing: in the presence 
of a strong laser pulse, the highly excited states acquire 
finite widths owing to ionization and the finite pulse du- 
ration, so that the electron moves in a quasi-continuum. 

In the case that we investigated (ag = 3, L = 0.3a), 
cf. Fig. ^ the enhancement as well as the general har- 
monic yield were markedly reduced compared with the 
untruncated potential. This is not so for a comparison 
of a realistic single-active-electron binding potential to a 
(regularized) zero-range potential in three spatial dimen- 
sions, neither for HHG nor for ATI 

In conclusion, we have confirmed the significance of 
resonant enhancements for high-order harmonic genera- 
tion. The resonances occur, when an appropriate highly 
excited state or, in the absence of such a state, the con- 
tinuum threshold are ponderomotively upshifted so that 
they become multiphoton resonant with the ground state. 
The mechanism is similar to high-order above-threshold 
ionization. In the latter case, remarkably, the multipho- 
ton resonance continues to be significant while the ion- 
ization process is already deeply in the tunneling regime 
|l9| . For high-order harmonic generation, this remains 
to be investigated. 
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TABLE I. Energies e„ of the field-free ground state and the first four excited states for the untruncated potential (g) as well 
as several truncated versions thereof, as considered in Fig. ^. No entry means that the state is no longer bound. Multiphoton 
resonances with the respective state occur for intensities rjn = integer + Arjn- 



untruncated ao = 8 L = 31.78 ao = 3 L = 31.78 ao = 3 L = 4.77 



n 


En 


Ar7„ 


£71 




En 




e„ 


Arjn 





0.7566 


0.04 


0.7566 


0.04 


0.7566 


0.04 


0.7565 


0.05 


1 


0.0846 


0.16 


0.0845 


0.16 


0.0825 


0.13 


0.0524 


0.74 


2 


0.0465 


0.67 


0.0453 


0.64 


0.0389 


0.55 






3 


0.0216 


0.33 


0.0124 


0.21 


0.0036 


0.09 






4 


0.0156 


0.25 


0.0014 


0.06 
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